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CALIBRATION OF [O IV] 26 AS A MEASURE OF INTRINSIC AGN LUMINOSITY 
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ABSTRACT 

We compare [O IV] 25.89 /im emission line luminosities with very hard (10-200 keV) X-rays from 
Swift, Integral, and BeppoSAX for a complete sample of 89 Seyferts from the Revised Shapley-Ames 
sample. Using Seyfert Is, we calibrate [O IV] as a measure of ACN intrinsic luminosity, for particular 
use in high-obscuration environments. With this calibration, we measure the average decrement in 
14-195 keV X-ray to [O IV] luminosity ratio for Seyfert 2s compared to type Is. We find a decrement 
of 3.1 ± 0.8 for Seyfert 2s, and a decrement of 5.0 ± 2.7 for known Compton-thick Seyfert 2s. These 
decrements imply column densities of approximately logNn = 24.6 cm^'^ and 24.7 cm~^, respectively. 
Thus, we infer that the average Seyfert 2 is more highly obscured and intrinsically more luminous 
than would be inferred even from the very hard X-rays. We demonstrate two applications of the 
hard X-ray to [O IV] ratio. We measure a column density for the extremely obscured NGC 1068 of 
XogNn = 25.3-25.4 cm~^. Finally, by comparing [O IV] luminosities to total infrared luminosities for 
twelve bright ultraluminous infrared galaxies, we find that four have substantial AGN contributions. 
Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — infrared: galaxies 



1. INTRODUCTION 

The optical hue [O III] 5007 A is widely used to 
measure the intrinsic luminosity of active galactic nu- 
clei (AGN) (iHeckman et al.|[2005[ ). However, there are 
cases, such as Seyfert 2s, ultra-luminous infrared galax- 
ies, and Compton-thick AGN, where the extinction along 
our line of sight to the narrow line regions may be high, 
and a more extinction -robust diagnostic may be needed. 
iMelendez et al.l ()2008D recently demonstrated a correla- 
tion in AGN between the [O IV] 25.89 /im emission line 
luminosity and the very hard X-ray luminosity (14-195 
keV) as measured by the BAT instrument on Swift. The 
sample of 40 AGN in Melendez et al., selected by hard 
X-ray-flux, was sufficient to demonstrate a correlation, 
but insufficient to calibrat e the relationship. In a com - 
panion paper to this one piamond-Stanic et al.ll2009D . 
we find that Seyfert Is and Seyfert 2s from a complete 
sample have statistically indistinguishable [O IV] lumi- 
nosity distributions, which further supports [O IV] as an 
AGN luminosity diagnostic. Here we compare the [O IV] 
and hard {E > 10 keV) X-ray luminosities of a complete 
sample of 89 Seyferts, bolometrically correct the hard X- 
rays to infer intrinsic AGN luminosities, and calibrate [O 
IV] as a measure of intrinsic AGN luminosity. 

One potential problem with [O IV] 26 fim is that its 
critical density of log ncrit = 4.06 cm~^ is lower than, for 
example, that of [Ne V] 14.32 /im (log n^rit — 4.70 cm~^), 
[Ne V] 24.31 /im (logric^t = 4.44 cm-^), or [O HI] 
5007 A ([ogHcrit = 5.8 cm~^)^. Aniano et al. (in prep.) 
find, using the same Seyfert sample as this paper, that 
the [Ne V] 14 / [Ne V] 24 fim line ratio is almost always 
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<2.0, which for temperatures of < 3 x 10^ K implies 
densities below 10"' cm^'^. Thus, [Ne V]/[Ne V] indicates 
that the densities are below the [O IV] 26 /im critical 
density, though this measurement applies only to the [Ne 
V]-emitting portion of the narrow-line region. 

Another potential problem, which affects any intrinsic 
luminosity diagnostic based on emission lines from the 
narrow-line region (NLR) , is that one must assume that 
the NLR enjoys an unobscured line of sight down to the 
central engine. If this assumption is violated, then the 
NLR will be ionized by less than the true continuum 
strength, and consequently NLR diagnostics will under- 
estimate the intrinsic AGN luminosity. 

For now, we develop [O IV] as a luminosity indicator, 
with the caveats that it may cease to work when densities 
exceed critical, and that like [O HI] 5007A it assumes 
the narrow line region has not been shielded from the 
full intensity of the central engine. In fjTl we examine 
[O IV] luminosities in ultra-luminous infrared galaxies 
(ULIRGs), which as deeply-obscured galaxies may be 
worst -case environments to test the [O IV] diagnostic. 

2. THE SAMPLE AND THE DATA 

Calibration of [O IV] against hard X-ray luminos- 
ity is best done within a complete sample to avoid 
flux biases. By contrast, the members of the Melen- 
dez et al. sample used to establish the correlation be- 
tween [O IV] and hard X-ray were selected as bright 
Swift detections, and therefore may well be biased to- 
ward X-ray-brightness. Our sample, which is limited 
by host galaxy magnitude, is the set o f 89 Se yferts from 
iMaiolino fc Rie kc (1995) and lHo et all ([1997D with Br < 
13. T his Revised Shapley-Ames (RSA. IShaplev &: AmesI 
19321) Seyfert sample is des cribed in more detail by 



Diamond-Stanic et al.l (120091). We take d i stance s and 
classifications from iDiamond-Stanic et al.l (|2009f) : the 
median distance is 22 Mpc. Optical spectral classifi- 
cations of Seyfert 1-1.5 we group as "Seyfert 1", and 
we group Seyfert 1.8-2 classifications as "Seyfert 2". 
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From column densities reported in the literature from 
X-ray measurements, we classify the AGN into reported 
Compton-thin, reported Compton-thick, and Seyferts 
whose column densities are not reported in the litera- 
ture. 

Spectra were obtained by the Infrared Spectro- 
graph (IRS) onboard Spitzer. We use [O IV] fluxes 
from the LLl order, as m easured and tabulated by 
iDiamond-Stanic et al.l ()2009[ ). These fluxes have been 
corrected for contamination fro m the [Fe II] 25.99 /im 
emissi on line, as described by IDiamond-Stanic et al.l 
(|2009( ). For the few AGN with undetected [O IV], we use 
the 3(7 upper limit. For the very bright sourc es NGC 1068 
and Circinus, we take ISO/SWS fluxes fromlSturm et all 
(l200l . and for NGC 4945 from lSpoon et al.l (|2000n . 

We take hard X-ray fluxes from published catalogs: 

• Swift B AT (14-195 keV) fl uxes from the 22 month 
survey (|Tueller et al.ll2009D : 



• BcppoSAX PDS (20-100 keV) fluxes from|Dadina 

(2007); and 



• Integral IBIS (17-60 

iKrivonos etHI (|2007f ). 



keV) fluxes from 



By using all three hard X-ray satellites, we obtain mean- 
ingful upper limits (since the Swift survey was all-sky), 
probe to fluxes fainter than the Swift all-sky detection 
limit, and gauge the importance of X-ray variability. 
Since Integral and BeppoSAX observations were tar- 
geted, the depth of coverage is highly non-uniform across 
the sky, so we do not estimate upper limits on non- 
detections. Since the Swift/BAT survey by design has 
fairly uniform all-sky coverage, for non-detections we plot 
upper limits of 3.1 x 10~^^ erg s~^ cm~^, which is the 
4.8cr depth the 22 month Swi ft survey reached for 90% 
of the sky ()Tueller et al.l[2009l ). 

We scale BeppoSAX fluxes to equivalent Swift fluxes 
using the bolometric corrections calculated in fj3]for the 
[Marconi et al. (2004) template; for AGN detected by 
both instruments the scaled fluxes agree within errors. 
We scale Integral fluxes by an empirical factor of xl.61, 
rather than the factor of xl.99 derived from the Mar- 
coni template. This empirical factor is the best-fit lin- 
ear relation for RSA Seyferts detected by both Integral 
and Swift. ^ The offset from the predicted relation is not 
surprising, since the Integral fluxes assume the spectral 
shape of the Crab nebula, which is softer than AGN over 
this spectral range (the Crab neb ula has F = 2.14 o ver 
the 1-700 keV range according to lKuiper et al.l[2001l .) 

3. THE BOLOMETRIC CORRECTION 

One wants to consider the total amount of "direct" 
or "intrinsic" AGN emission, which is radiated primarily 
in the X-ray, UV, and optical. This is the luminosity 
inferred to have been emitted into An steradians before 
any obscuration. One should separately consider emis- 
sion that has been absorbed and re-radiated into the in- 
frared, since this will depend on the geometry of the 
absorbing material. A so-called bolometric correction, 
based on model spectral energy distributions, is used to 
translate a luminosity measured at a particular band (for 

® After dropping the two weakest detections. 



example, hard X-rays) into the intrinsic luminosity of the 
AGN. (This is termed a "bolometric correction" in the 
literature, even when, as in this case, the reprocess ed in- 
frared emission is deliberately ignored.) Marconi et al.l 
(|2004!) constructed a template spectrum using broken 
power laws in the optical-UV, a Rayleigh-Jeans black- 
body tail in the near-IR, and a power-law spectrum with 
F = 1.9 plus a reflection component in the X-ray. The 
X-ray luminosity was scaled t o match the luminosi ty de- 
pendence of aox measured bv lVignali et al.l (|2003f ). The 
resulting template s have lowe r X-ra y-to-intrinsic lumi- 
nosi ty ratios than lElyis et al.l (|1994f ). as expected since 
the I Elvis et all ()1994[ ) sample was c hosen to be X-ra 



bright. In correcting for that bias, lElvis et "aTl ()200: 
shifted the X-ray portion of their template down 0.135 
dex; this accounts for all bu t 0.0 5 dex of the off s et be - 
tween the lElvis et all (|1994D and IMarconi etHI (|2004D 
templates. 

From these templates, Marconi et al. calculate the con- 
version from L(2-10 keV) to L(intrinsic) as a function 
of L(intrinsic), and plot the relation in their figure 3b. 
We simply scale this conversion to higher X-ray energies. 
Fitting a fifth-order polynomial to the Marconi et al. 
template and integrating over the relevant X-ray bands, 
the scaling relations are: 

/(14- 195 fcey)//(2- 10 fcey) = 2.67 (1) 

/(20 - 100 keV)/,f{2 - 10 keV) = 1.74 (2) 

/(17-60 fcey)//(2-10 fceF) = 1.34 (3) 

for the Swift BAT, BeppoSAX PDS, and Integral IBIS 
bands, respectively. We then fit the dependence of in- 
trinsic luminosity L/ on 2-10 keV luminosity L2-10 in 
the Marconi models as: 

log(L//Lo) = 0.03776 [log(i2-io/io)]' + (4) 
0.5340 log(L2-.io/£o) + 2.276 

This equation, along with the scaling relations in eqn. 
1-3, can be used to infer an intrinsic AGN luminosity 
from a measured high - energ y X-ray luminosity. 

The IMarconi et al.l (j2004l ) template is well-suited to 
this application because it ac counts for the know n lu- 
minosity dependence of aox (jVignah et al.l 120031 ) . and 
because it e xtends in energy up to 1 MeV, unlike for 
example the lElvis et al.l (|1994f ) template which stops at 
40 keV . We have already shown that the IMarconi et al.l 
(|2004f l template agrees well with the revised Elvis et al. 
template (2002), and will thus generate similar bolomet- 
ric correctio ns. We now check t h e bol ometric correc- 
tion against iVasudevan fc FabianI (12007). W e consider 
the 54 AGN from IVasudevan fc FabianI (|2007D that have 
intrinsic AGN luminosit ies estimated from FUSE. Six- 
teen are detected in the iTueller et al.l (|2009( ) Swift cat- 
alog. In Figure [T] we plot the AGN luminosity from 
IVasudevan fc FabianI (|2007f) against the Swift BAT lu- 
minosity. We plot a simple linear fit considering only 
detections, as well as a non-linear fit that includes BAT 
upper limits using the Buckley-James regression method 
(dsobc et al. 1986). Fi gure [1] shows that the non-linear 
fit agrees extremely well with the Marconi relation. 

4. COMPARING [O IV] AND HARD X-RAY LUMINOSITIES 
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Fig. 1. — Comparison of bolometric luminosity (inferred from 
FUSE f ar-UV) and Swift BAT hard X-ray luminosity for the 16 
AGN in IVasudevan fc FabianI H2007I ) detected by BAT. Two fits 
are plotted: a proportional fit to the detections (dashed line) of 
-L//Lq= 4.42 Lbat; a-nd a non-linear fit that includes the BAT 
upper limits using the Buckley- James regression method ( grey solid 
line) o{log{Li / Lq) = 1.3A5log{LBAT/LQ) — 2.60. The non-linear 
fit ag rees extremely well with the relation from IMarconi et al.l 
P004l ^ (black solid line), where a bolometric correction has been 
applied between BAT and 2-10 keV. 

Figure [2] plots Swift-equivalent luminosity against [O 
IV] luminosity. Each detection is plotted separately, with 
symbols to indicate the satellite (squares for Swift, cir- 
cles for BeppoSAX, triangles for Integral). (For example, 
an AGN detected by all three satellites is represented by 
a cluster of three points at the same L[0 IV] value.) The 
AGN are split into four panels using classifications de- 
scribed in ij2l Seyfert 1; Compton-thin Seyfert 2; Seyfert 
2 with unknown column density; and Compton-thick 
Seyfert 2. 

Next, we use the bolometric corrections of 
IMarconi et al.l (|2004f ) to convert the Swift-equivalent 
X-ray luminosities into intrinsic AGN luminosity, and 
plot against L[0 IV] in Figure [H Symbols indicate 
column density and spectral classification. We confirm 
the result of Melendez et al. (2008) that Seyfert 2s have 
systematically lower X-ray/ [OIV] luminosity ratios than 
Seyfert Is. In addition, with our larger sampl e of 18 
Comp ton-thick AGN (compared to 4 in Melende z et al.l 
UdoI)), of which 15 were detected in [O IV] , and 6 were 
detected in hard X-ray, we show that Compton-thick 
AGN have still lower ratios. 

For each AGN class, we fit the ratio of intrinsic AGN 
luminosity inferred from the hard X-rays to [O IV] lumi- 
nosity using t he Kaplan -Meier estimator of the survival 
function (Fei eelson fc Nelson 1985), as imp lemented in 
the ASURV package p^Vallev et al.lll992( ). This esti- 
mator incorporates both detections and upper limits to 
provide a maximum likelihood reconstruction of the true 
distribution function, even in cases where there are more 
upper limits than detections (as for the Compton-thick 



AGN.) 

For AGN detected by multiple X-ray instruments, we 
use the average scaled X-ray flux. For Seyfert Is, the 
best-fit ratio is 2550 ± 740 (with a scatter of 0.4 dex); 
for Seyfert 2s it is 810 ± 170; for known Compton-thin 
Seyfert 2s it is 1320 ±300; and for known Compton-thick 
Seyfert 2s it is 510 ±270.^ If we exclude the anomalously 
X-ray-bright NGC 4945 from the fit, the best-fit ratio for 
the Compton-thick Seyfert 2s falls to 240 ± 70.^ Thus, 
Seyfert 2s have an average hard X-ray to [O IV] ratio 
that is 3.1 ± 0.8 times lower than Seyfert Is. Known 
Compton-thin Seyfert 2s show a decrement of 1.9 ± 0.5, 
and known Compton-thick Seyfert 2s show a decrement 
of 5.0 ± 2.7 (excluding NGC 1068.) 

The two-sample ASURV statistical tests find a < 
0.06% probability that the X/[0 IV] ratios for Seyfert 
Is and Seyfert 2s are drawn from the same population. 
It is similarly improbable (< 0.09%) that Seyfert Is and 
Compton-thick Seyfert 2s are drawn from the same pop- 
ulation. However, the offset between Sis and Compton- 
thin Seyfert 2s is of lower statistical significance (9% 
probability of being drawn from the same population) . 

5. ESTIMATING COLUMN DENSITIES 

We now interpret the difference in hard X-ray/ [O IV] 
ratio between type 1 and type 2 Seyferts. One expla- 
nation would be that the [O IV] luminosities for the 
type 1 Seyferts are biased high. In a companion paper 
([Diamond-Stanic et all (2009) , we find no significant dif- 
ferences between the [O IV] luminosity distributions of 
type 1 and type 2 Seyferts; the same conclusion can be 
reached by inspecting Figures O and [3] Thus, there is no 
evidence for a bias toward high [O IV] luminosities. 

Therefore, we use L[0 IV] as a measure of intrinsic 
AGN luminosity, and interpret the offset in hard X- 
ray/[0 IV] ratios as caused by substantial absorption 
in the E>10 keV spec tra of Seyfert 2s, as suggested by 
iMelendez et~all (|2008[ ). At such high energies, the cross- 
section for Compton scattering is larger than that of pho- 
toelectric absorption. Consequently, one expects a typi- 
cal high-energy photon to experience multiple Compton 
scatterings, losing energy each time to electron recoil, 
until the photon either escapes, or loses sufficient en- 
ergy and traverses a sufficient path that it is absorbed 
photoelectrically. This makes the problem relativistic, 
optically-thick, and geometry-dependent, and as such, 
Monte Carlo simulations are required to i nterp ret the 
flux decrement. Fortunately, iMatt et al] (|1999( ) simu- 
lated exactly this situation. We numerically integrate 
their spectra over the Swift/BAT, BeppoSAX/PDS, and 
Integral/IBIS energy bands and compare to the injected 
spectrum; the resulting flux decrements are listed in Ta- 
blelB 

In the IMatt et al.l jTMl) model, suppressing a BAT 
flux by a factor of 3.1 ±0.8 (as observed for all Seyfert 2s 
in our sample) requires a column density of log N(H) = 
24.6lo;2 cm~^. Suppression by a factor of 5.0 ± 2.7 (as 
observed for known Compton-thick Seyfert 2s®) requires 
log N{H) = 2A.7t°oi cm . These results are consistent 

^ We exclude NGC 1068 from the fit because its ratio is so low, 
and consider it separately in ^ 
* Excluding NGC 1068 as well. 
9 Excluding NGC 1068. 
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wi th the column density o f log N{H) = 24.5±0.1 inferred 
by iMelendez et al.l ()2008D for their sample of 17 Seyfert 
2s. Thus, the [O IV] to hard X-ray ratio confirms that, 
on average, Seyfert 2s are more obscured in the hard X- 
rays than Seyfert Is, and predicts the obscuring column 
densities expected for Compton-thick AGN. 

Surprisingly, the average ratio for the known Compton- 
thin Seyfert 2s also suggests high columns. Suppression 
by a factor of 1.9 ± 0.5 can be explained by log N{H) = 
24.3^03 cm~^. By contrast, the median column den- 
sity reported in the literature for these AGN, based 
on 2-10 keV measurements, is log N{H) ~ 23.0 cm^^ 
(jPiamond-Stanic et al.l 120091 ). As such, the hard X- 
ray/[0 IV] ratios imply higher columns (Compton-thick 
or nearly so) than inferred from the 2-10 keV spectra. 
However, this is not the only plausible interpretation. 
As discussed in ^ the offset in hard X-ray/ [O IV] ra- 
tio between Compton-thin Seyfert 2s and Seyfert Is is 
of low statistical significance in our sample. As such, we 
cannot rule out Compton-thin Seyfert 2s having, on av- 
erage, low columns. Table [1] shows that column densities 
below log N{H) ^ 24 cm^^ should produce only a small 
flux decrement. Accordingly, for such low columns, soft 
X-ray spectra should return a more precise measurement 
of the column density. 

Figure [3] does show that a substantial minority of the 
Compton-thin Seyfert 2s have X/[0 IV] ratios as low as 
known Compton-thick AGN. In order of increasing X/[0 
IV] ratio, these are NGC 1365, NGC 7582, NGC 2992, 
NGC 7314, NGC 3081, and NGC 5728. i° There several 
possible explanations for these low ratios: they may have 
nearly Compton-thick or Compton-thick columns, have 
column densities that vary dramatically with time, or 
have dramatic variations in hard X-ray flux. The litera- 
ture reveals that all three explanations appear to be at 
work: 

• NGC 1365 varies between Comp ton-thin and 
Comp ton-thick in a matter of days (|Risaliti et al.l 
[2001: 

• NGC 7582 has experienced dramatic changes in its 
2-10 keV flux c onsistent with a variab le Compton- 
thick absorber (jPiconceUi et al.ll2007b : 

• NGC 2992 was measured to have a column of only 
7x lO^^cm-^ (Colbert et al. 2005), but its 2-10 keV 
flu x has varied by a factor of 20 over twenty years 
(GiUi et al.|[2000i). and its BAT flux has varied by 
a factor of 6 ( Beckmann et al.ll2007f ): 



• NG C 7314 was me asured to have a low column den- 
sity (jRisalitil 12001 ): 

• NGC 3081 was m easured to have a colu mn density 
of 6 X 10^3 cm-2 UMaiolino et al.lll998D. and varies 
in th e BAT band by a factor of 7 ( Beckmann et al.l 
[2003); 

• NGC 5728 is borderline Compton thick with a col- 
umn of 8 X 10^3 cm-2 (Zhang et al. 200^. 

The remaining Compton-thin Seyfert 2s have an average hard 
X-ray to [O IV] offset from the Seyfert Is that is of even lower 
statistical significance. 



Thus, five of these six AGN have a plausible explanation 
for their low X/[0 IV] ratios; in four cases the explana- 
tion appears to be that the column is almost Compton- 
thick, or transitions between Compton thick and thin. 
This comparison also suggests that an inherent drawback 
of calibrating [O IV] against the hard X-rays is that the 
hard X-ray flux may vary with time. Still, the system- 
atically low X/[0 IV] ratios observed for all Seyfert 2s 
in this sample strongly suggest that the typical Seyfert 
2 has a signiflcantly attenuated hard X-ray flux, and is 
considerably more powerful than would be inferred from 
the hard X-rays. 

6. APPLICATION: NGC 1068 

The famous Compton-thick NGC 1068 is the most 
highly-obscured AGN in our sample, based on its hard 
X-ray to [O IV] luminosity ratio. Its [O IV] luminosity^^ 
of 1.2 X 10* L© predicts an intrinsic AGN luminosity 
of (3.1 ± 0.9) xlO^-^ Lq using the Seyfert 1 cahbration. 
The JO IV] flux is take n from high resolution ISO spec- 
tra (jSturm et al.ll200l . so [Fe II] contamination is not 
an issue; in any case the [Fe II] line strength is only 4% 
that of [O IV]. The hard X-ray to [O IV] ratio is 490 x 
lower than the Seyfert 1 best flt, implying a c olumn den- 
sity o f log TV// = 25.3 to 25.4 cm'^ in the llVEatt et all 
()1999l ) models. (The column should be still higher if 
scattered light contributes substantially to the observed 
X-ray flux.) This inferred column is consistent with the 
lower limits inferred through a similar comparison to [O 
III] 5007A ()Matt et aLllTglTl [2001 . In a future paper, 
we plan to compare this inferred intrinsic luminosity in 
detail with the reprocessed mid- infrared SED. 



7. 



APPLICATION: ULTRA-LUMINOUS INFRARED 
GALAXIES 



As an application of the [O IV] calibration, we con- 
sider the [O IV] luminosities for twelve nearby ultra- 
luminous inf r ared gala xies (ULIRGs ) , usin g data from 
lArmus et all (|2004l) and lArmus et al.l (|2007D . which com- 
prises the ten brightest ULIRGs plus two others. Only 
five ULIRGs have detected [O IV]: IRAS 05189-2524, 
IRAS 13428-f5608 (Mrk 273), IRAS 09320-1-6134 (UGG 
5101), Mrk 1014, and Mrk 463e. Using the relation es- 
tablished above for Seyfert Is, we infer intrinsic AGN lu- 
minosities from L[0 IV] and compare to the total infrared 
(8-1000 Aim) luminosities. The inferred L(AGN)/L(IR) 
percentages are 49%, 97%, 22%, 170%, and 520%, re- 
spectively. (A 30% errorbar should be applied for the 
uncertainty in the Seyfert 1 calibration.) For ULIRGs 
with undetected [O IV] we infer percentages below 8%. 

Based on excitation diagrams using [Ne V]/[Ne II], 
[O IV]/[Ne II] , and the equivalen t width of the 6.2/im 
PAH feature, lArmus et"aI1 (|2007f) picked out three of 
these (Mrk 463e, Mrk 1014, and IRAS 05189-2524) as 
being likely AGN-dominated. In addition, they esti- 
mate Mrk 273 has an AGN contribution of 50-70%. 
UGG 5101, which [O IV] indicates has a modest AGN 
con tribution, does not stand out as AGN-dominated in 
the lArmus et al.l (|2007l) excitation diagrams; however, it 
does ha ve detected Fe K a w ith an equivalent width of 
400 eV (jlmanishi et al.|[2003[ ). which signals that an ob- 
scured AGN is present and may possibly be bolometri- 

Assuming a distance of 14.4 Mpc. 
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cally important. Thus, the L[0 IV]/L(IR) ratio gives re- 
suhs consistent with the excitation diagrams as to which 
ULIRGs are AGN-dominated. 

Does the normahzation of L[0 IV] to L(AGN), cal- 
ibrated from Seyfert Is, appear roughly correct for 
ULIRGs? The preliminary answer appears to be yes, 
since the inferred AGN fractions are of order unity for 
four of the five the [O IV] -detected ULIR Gs, consistent 
with t he results of the excitation diagrams. lArmus et al.l 
(|2007D did note that the AGN fractions implied by 
[Ne V]/[Nc II] or [O IV] /[Ne II] line ratios were sig- 
nificantly lower than implied by the mid-IR slope or 
the 6.2 /xm PAH feature equivalent width; we do not 
see this effect for L[0 IV]/L(IR). The one clear outlier 
in L[0 IV]/L(IR) is Mrk 463e, for which we infer an 
L(AGN)/L(IR) ratio significantly above unity. Such a 
result is physically plausible if most of the AGN energy 
is not reprocessed into the infrared, but escapes via the 
X-ray/UV/optical. Indeed, the infrared luminosity for 
Mrk 463e is only 5 x 10" Lq, below the ULIRG cutoff, 
but its total luminosity is considerably hig her, well above 
the ULIRG threshold (" Armus et all 12004). In addition, 
Mrk 463e has remarkably weak PAH features, which ar- 
gues that the starburst contribution is small. Thus, given 
that Mrk 463e is apparently an AGN-dominated ULIRG 
with a low ratio of reprocessed to intrinsic light, the high 
measured ratio of L[0 IV]/L(IR) seems reasonable and 
physical. 

It remains unclear how to interpret the [O IV] non- 
detected ULIRGs, especially those like Mrk 231, which 
is a broad-line AGN. One explanation is that these galax- 
ies are indeed dominated bolometrically by star forma- 
tion, not AGN. An alternate explanation would be that 
they have bolomctrically-important AGN, but that the 
[O IV] diagnostic has failed, perhaps because densities 
exceed the [O IV] critical density, or because the nar- 
row line regions have been shielded from the central en- 
gine. Thus, while the [O IV] diagnostic appears to work 
well for Seyferts in the RSA sample, and identifies what 
are thought to be the most AGN-dominated of bright 
ULIRG, more work is still required to understand its 
strengths and limitations of the diagnostic in extremely 
obscured environments. 

8. CONCLUSIONS 

We present the first calibration of [O IV] 25.89 fim as 
a measure of intrinsic AGN luminosity. The [O IV] line 
is bright, is detected in hundreds of AGN in the Spitzer 
archive, and will be detectable to Herschel and JWST. 

We find that Seyfert 2s show systematically lower 
{E > 10 keV) ratios of X-ray to [O IV] luminosity than 
Seyfert Is, by a factor of 3.1 ± 0.8. For AGN previ- 
ously identified as Compton-thick, the observed decre- 
ment is 5.0 ± 2.7. We interpret this as caused by ab- 
sorption that affec ts even the very h ard X-rays. Monte 
Carlo simulations (jMatt et al.lll999f) associate these X- 
ray/[0 IV] flux decrements with average column densi- 
ties of logA^(iJ) = 24.6tg:2 and 24.7t^:3 cm"^ re- 
spectively. This, we find substantial obscuration in most 
Seyfert 2s. 

We briefly explore use of [O IV] 26 /im to estimate 
AGN power even in very obscured systems such as 
Compton-thick AGN and ULIRGs. For Compton-thick 
AGN we infer sensible average column densities, and 



are able to infer a column density even for the ex- 
tremely obscured NGG 1068. Application to ULIRGs 
is less straightfoward, since only five of the twelve bright 
ULIRGs we examine have detected [O IV]. For [O IV]- 
detected ULIRGs, the L(0 IV) /L(IR) ratios imply signif- 
icant AGN contribution. The question remains whether 
the non-detections have lower AGN fractions than the de- 
tections; it is possible that the [O IV] diagnostic breaks 
down in the extreme environments of some ULIRGs. 

We conclude by noting the recent interest in using very 
hard X-rays to find obscured AGN, due to their pene- 
trating power. Several sensitive, high-resolution £^ > 10 
keV wide-field or pencil-beam X-ray satellites are now 
in preparation or planning, with scientific goals of find- 
ing even the most obscured AGN in the nearby universe, 
and identifying the population of Compton-thick AGN 
presumed to contribute to the hardness of the X-ray 
background. Our results suggest that even the 20-200 
keV band suffers considerable flux attenuation in typi- 
cal Seyfert 2s, especially those that are Compton-thick. 
Thus, our results predict that X-ray surveys will probe 
considerably smaller volumes for obscured AGN than for 
unobscurcd AGN, even when very hard energy bands are 
used. The [O IV] 26/im line may play a critical role in 
measuring this absorption, and determining the true lu- 
minosities of even highly-obscured AGN. 

Acknowledgments: We thank G. Matt for kindly pro- 
viding Monte Carlo simulation results from iMatt etldl 
(1999;) for additional column densities. JRR was sup- 
ported by a Spitzer Space Telescope Postdoctoral Fel- 
lowship. 
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TABLE 1 

Ratio of em ergent to input ha rd X-ray flux, from 

IMATT et AL.I (119991) MODELS. 



log Nh (cm-2) Swift/BAT SAX/PDS Integral/IBIS 



23.0 


0.96 


0.96 


0.96 


23.5 


0.88 


0.89 


0.87 


23.8 


0.78 


0.80 


0.77 


24.0 


0.69 


0.71 


0.68 


24.1 


0.64 


0.66 


0.63 


24.2 


0.58 


0.61 


0.57 


24.3 


0.51 


0.54 


0.50 


24.4 


0.44 


0.48 


0.44 


24.5 


0.37 


0.41 


0.38 


24.6 


0.30 


0.34 


0.31 


24.7 


0.23 


0.27 


0.25 


24.8 


0.17 


0.20 


0.19 


24.9 


0.11 


0.13 


0.13 


25.0 


0.063 


0.077 


0.074 


25.1 


0.048 


0.058 


0.056 


25.2 


0.028 


0.034 


0.033 


25.3 


0.0032 


0.0039 


0.0039 


25.4 


0.0016 


0.0020 


0.0020 


25.5 


0.00010 


0.00013 


0.00013 
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Fig. 2. — Hard X-ray luminosities as a function of [O IV] line luminosity. Seyfert 1 AGN are plotted in the top panel; known Coinpton-thin 
Seyfert 2 in the seeond panel; Seyfert 2 of unknown column density in the third panel; and Compton-thick Seyfert 2 in the bottom panel. 
Swift BAT detections arc squares (upper limits are black partially— filled arrows); BeppoSAX PDS detections are circles (upper limits are 
black filled arrows); and Integral detections are triangles (upper limits are grey arrows.) BeppoSAX and Integral fluxes have been scaled 
to equivalent fluxes in the Swift BAT bandpass. 
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Fig. 3. — Intrinsic AGN luminosity, as inferred from the hard X-rays, versus [O IV] luminosity. Seyfert 1-1.5 AGN are plotted in black 
as squares (detections) and arrows (non-detections); Compton-thin Seyfert 1.8—2 AGN are plotted as red open circles and arrows; Seyfert 
1.8—2 of unknown column are plotted as blue arrows or squares; and Compton-thick AGN are plotted as light grey open squares and arrows. 
The best fit X-ray-inferred intrinsic to [O IV] luminosity ratio is plotted for each AGN class, in the same color scheme; the fits are given 
in m 



